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Abstract 

We study the impact of a strongly first-order electro-weak phase transition on the thermal relic 
abundance of particle species that could constitute the dark matter and that decoupled before 
the phase transition occurred. We define a dilution factor induced by generic first-order phase 
transitions, and we explore the parameter space of the minimal supersymmetric extension to the 
Standard Model to determine which phase transition temperatures and dilution factors are rele- 
vant for the lightest neutralino as a dark matter candidate. We then focus on a specific toy-model 
setup that could give rise to a strongly first-order electro-weak phase transition, and proceed to 
a detailed calculation of dilution factors and transition temperatures, comparing our findings to 
actual neutralino dark matter models. Typical models that would produce an excessive thermal 
relic density and that can be salvaged postulating a strongly first-order electro-weak phase transi- 
tion include massive (multi-TeV) wino or higgsino-like neutralinos, as well as bino-like neutralinos 
in a wider mass range, with masses as low as 400 GeV. If LHC data indicate an inferred ther- 
mal neutralino relic abundance larger than the cold dark matter density, the mismatch could thus 
potentially be explained by electro-weak scale physics that will also be thoroughly explored with 
collider experiments in the near future. 
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I. INTRODUCTION 



The particle content of the Minimal Supersymmetric extension of the Standard Model 
(MSSM) of particle physics (see e.g. [1]) provides a stable weakly interacting massive particle 
which is, in principle, a viable particle dark matter candidate: the lightest neutralino. If 
the latter is the lightest i?-parity odd particle, its thermal relic abundance can be close to 
the inferred density of dark matter on cosmological scales whose ratio to the universal 
critical density is estimated to be f^astro = Pdm/Pc = 0.113 [3], with h the present day 
Hubble expansion rate normalized in units of 100 km per sec per Mpc. 

Several studies have pointed out that the thermal relic abundance of neutralinos in MSSM 
models where they are the lightest supersymmetric particles actually ranges over several 
orders of magnitude, simple estimates of their relic abundance being therefore only a generic 
order-of-magnitude estimate (see e.g. |4, Isl). Specifically, if the lightest neutralino has 
unsuppressed couplings to gauge bosons, and is heavy enough for annihilation to massive 
weak interaction gauge bosons to be kinematically open, it rapidly annihilates into VI^+VI^~ 
and/or ZZ final states, and its relic density f^particie is below the dark matter density, at least 
for neutralinos lighter than about a TeV. If the lightest neutralino has suppressed couplings 
to gauge bosons, as is the case for bino-like or singlino-like neutralinos, or if neutralinos are 
very heavy, the neutralino relic density is instead typically much larger than the universal 
dark matter density. 

MSSM models with under- abundant neutralino relic density are phenomenologically per- 
fectly viable: the lightest neutralino can very well not be the only contributor (or it can 
be a sub-dominant contributor) to the universal dark matter, in a standard cosmological 
setup. While modified cosmological setups can be concocted to enhance the thermal relic 
abundance of neutralinos in those cases (for instance in the presence of non-thermal pro- 
duction 15,1, or of a faster expansion rate at the epoch of the lightest neutralino freeze-out 

riQ 

|6|, lZI]), MSSM models with over-abundant relic neutralinos are, in principle, ruled out. A 
caveat to this conclusion is the possibility of an episode of entropy injection occurring at 
temperatures below the decoupling of neutralinos from the universe's thermal bath (entropy 
injection before freeze-out would not change today's relic abundance of neutralinos). 

Ref. [5] gives a thorough discussion of such a class of scenarios for the specific case 
of an additional scalar field that drives both the entropy dilution and possibly the non- 
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thermal production of neutralinos, as well as a modified cosmological expansion rate if the 
energy density associated to the scalar field dominates the universal energy budget. Specific 
examples of setups that can produce a dilution in the relic abundance of neutralinos (or, 
for that matter, of any other thermal relics decoupled from the thermal bath before entropy 



injection) include for example models with moduli 
with low-temperature or even weak-scale infiation [l 



11 



or Q-ball decays jol], and scenarios 



121]. While a phase transition 



that plausib! 



y occurred after the neutralino freeze-out is the QCD phase transition, lattice 



simulations [l3|] clearly indicate that the transition is not strongly first-order (see e.g. the 



discussion in Ref. [14]). 

In the present study, we consider a specific instance of a possible source of dilution for 
the thermal relic abundance of relic neutralinos from the early universe, namely the electro- 
weak phase transition (EWPT). In the minimal formulation of the Standard Model Higgs 
sector, the EWPT is very weakly first-order, or entirely absent, as found e.g. in the non- 
perturbative analysis of Ref. [15]. However, a strongly first-order phase transition can occur 
within the MSSM, or even in very simple extensions of the scalar sector of the Standard 

n 

Model (for instance via the inclusion of additional singlet scalars PJl]). Other possibilities 
include for instance models with multiple hidden sector scalars coupled only to the Higgs 
sector (which thus acts as a "portal" to the hidden sector, as envisioned in Ref. [iT, 18, 19]). 
This possibility, which includes scenarios with tree-level conformal invariance where the 
Higgs mass is generated via dimensional transmutation, has been discussed and studied 
in detail in Ref. [20] and [2l|]. A strongly first-order phase transition needs to actually 
be posited in the context of scenarios where the baryon asymmetry is produced at the 
EWPT, to prevent the "washout" of the generated baryon number density by sphaleron 
processes (for a pedagogical introduction to the framework of electro- weak baryogenesis see 



e.g. Ref. [221]; classic studies on MSSM baryogenesis include [23] and [24] ; recent discussions 
of the phenomenology of electro- weak baryogenesis in the MSSM are given in Ref. [251. 1261]). 
In this case, the order parameter can be considered to be the ratio of the SU(2) Higgs 
background field (j) = ^/2{H^) to the critical temperature for the phase transition, Tc, and a 
strongly first-order phase transition corresponds to ^c/Tc ^ 1 j22| . 

The entropy injection produced in any first-order phase transition can play a relevant role 
in the thermal history of species that froze out prior to the phase transition, since their relic 
density will be diluted away by an amount dependent upon the relative entropy injected to 
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the entropy in the species in thermal equihbrium. Interestingly, the CERN Large Hadron 
Collider (LHC) might potentially inform us on both the details of the scalar electro-weak 
sector and of the EWPT in particular (including possible non-minimal extensions of the 
scalar sector [igI), as well as on the mass of a putative dark matter particle. In some cases, 
data from the LHC might even enable us to infer, via the knowledge of particle masses and 
couplings, the thermal relic abundance of a stable neutralino that might be the main dark 
matter constitue nt ( explicit examples with realistic assumptions on the LHC performance 



are given in Ref. 



271]). A mismatch between the observed cold dark matter abundance and 



the inferred relic abundance of candidate particles discovered with colliders can have very 
profound consequences. If the resulting inferred thermal abundance were larger than the 
cold dark matter density, then the answer to the ensuing conundrum might lie in the EWPT 
and in the Higgs sector, and thus again in physics that can, and will, be tested with the 
LHC. 

The present study is organized as follows: in the next section we describe in detail how 
a first-order phase transition affects the abundance of a species that froze out prior to 
the temperature at which the phase transition occurred; the following sec. [Till explores the 
ranges of freeze-out temperatures and of relic abundances relevant to the case of the lightest 
neutralinos of the MSSM. The ensuing sec. HVl describes simple models for the EWPT, and 
sec. |V] describes their impact upon the relic density of MSSM neutralinos. These last two 
sections present our main results. Finally, sec. I VII summarizes and concludes. 



II. RELIC DENSITY DILUTION: THERMODYNAMICS 

We are concerned here with determining the dilution of a relic species that decoupled (or 
''froze-out" ) from the thermal bath in the early universe by entropy injected during a first- 
order phase transition occurring after the species' freeze-out. The basic thermodynamics of 
the dilution, which we describe in the present section, is fully general, and does not depend 
upon the specifics of the model for the phase transition. 

Let ^(0, T) be the finite temperature effective potential of the early universe, where (j) is 
the order parameter of the phase transition, for instance the vacuum expectation value of 
the neutral component of the Higgs field. We assume that, for the temperatures under con- 
sideration here, there are at most two minima of the potential. At very high temperatures. 
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the potential only has one minimum, at = 0. As the temperature drops, a second mini- 
mum develops with a corresponding effective potential value larger than that at = 0. The 
value of the potential at the second minimum relative to the first decreases with decreasing 
temperature so that at T = the second minimum is the absolute minimum and there is 
spontaneous symmetry breaking. We define a "critical temperature" as the temperature 
T = Tc when the two minima are degenerate. At some temperature < T^^ the system 
transitions from = to the new minimum: the tunneling probability from the = 
minimum to the true vacuum is of order unity. 

First, let us assume that the transition temperature is very close to the critical tem- 
perature Tc so that supercooling is negligible. Here, we follow the discussion of Megevand 
and Sanchez (see also j29|). Let 5+ be the entropy density of the high-temperature 
phase, and let 5_ be the entropy density of the low-temperature phase. We may then write 
the total entropy density as 

s = s+- fAs, (1) 

where As = 5+ — 5_ and / is the volumetric fraction of the system in the low temperature 
phase. Since the minima are degenerate at Tc, the system is in equilibrium and the total 
entropy is conserved. The entropy density then scales as 

. = ..(^;, (2) 

where is the scale factor of the universe at the beginning of the phase transition. The 
transition is complete once / = 1. Plugging this into Eq. ([T]) and combining with Eq. ([2]), 
one finds that the total expansion during the phase transition is 



ai J 1 — As/s^ 5_ 

This is the equation we use to find the dilution in sec. lIVBi 

Realistically, the system would not be in exact equilibrium — friction and collisions in the 
walls of bubbles of true vacuum would release entropy, slightly increasing the dilution js^ . 
This can be quantitatively captured by assessing the variation of the dilution factor with 
a variation to the temperature at which the phase transition occurs, this variation being 
driven by the mentioned effects (bubble collisions, friction in the bubble walls etc.). If the 
entire transition occurs at a temperature an amount AT below Tc, then one can show that 
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the change in dilution D is AD = —3D{D — 1)(AT/Tc), assuming a radiation dominated 
energy density. 

If "^Tc^ then the two phases are not in equihbrium at the beginning of transition and 
Eq. ([2l) does not hold. Instead, there are three distinct stages to the phase transition: a 
supercooling stage, a reheating stage, and a phase-coexistence stage. During supercooling, 
entropy is conserved, so that 

where is the scale factor at the minimum of supercooling. Assuming that reheating 
happens quickly relative to the expansion rate, the energy density p of the universe does not 
change during reheating. The entropy, however, does. If a large enough amount of reheating 
occurs, the universe will reach a phase coexistence stage at T = Tc. Conservation of energy 
then gives the initial fraction of the universe in the low-temperature phase at the beginning 
of phase coexistence: 

p+in) = p+(r,) - /o[p+(r,) - p_(r,)] 

(5) 

= P+{Tc) - foL 

JO — , loj 

where L is the latent heat of the transition at Tc. The entropy density during phase- 
coexistence is then 

.= (.+ -/oA.)(^)'. (7) 
Combining this with Eq. ([T]), we have the expansion during phase coexistence 



a/V_ l-foAs/s+{Tc) 
l-As/s+{Tc) ' 



which gives a total expansion of 



«/^'_ fl-foAs/s+{Tc)\ s+(T,) 



aij V 1- As/s+(re) ; s+{T,y 

We use this equation in sec. IIV CI where the transition is strongly first-order. 



(8) 



(9) 



A. Determining the transition temperature 



We specify here the definition we adopt for the transition temperature T*, given its 
relevance in determining the dilution factor for a given effective potential. Shortly after 
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the universe cools below the critical temperature, bubbles of true vacuum nucleate via 
thermal tunneling. Most of these are too small to grow — the pressure difference Ap = —AT 
between the true and false vacuum is not large enough to overcome the surface tension of 
their walls, so they collapse. Only large bubbles can grow. As the universe further cools, 
the nucleation rate of larger bubbles increases dramatically. The phase transition begins 
once the probability to nucleate a supercritical bubble in one Hubble volume is of order 1. 
Tunneling in cosmological phase transition was originally discussed in the seminal work of 
Ref. 3l|, l32| (see also Ref. 291]). For a pedagogical review of cosmological phase transition, 
see e.g. Kolb and Turner [331]. 

The tunneling probability per unit time per unit volume goes with temperature roughly 



as r ^ T^exp ^^/^^ where Ss is the three-dimensional Euclidean action 

2 



S. = 47r 



r dr 



(10) 



and where we assume spherical symmetry. The bubble shape (f){r) comes from the corre- 
sponding Euclidean equation of motion 

d 



2 # 
dr'^ r dr 



(11) 

0. Finally, the requirement 



and ^ „ 

dr \r=0 



with the boundary conditions lim^^oo 0(^) 
that one supercritical bubble nucleates per horizon volume in a Hubble time yields a phase 
transition temperature such that 



^ dT f2CMpi 



exp 



-S3(T)/T 



0{1) 



(12) 



where Mpi is the Planck mass, ( = and g is the effective number of relativistic 

degrees of freedom (see Ref. [22]). For temperatures around the electroweak scale, Eq. ( |T2l ) 
implies that S3/T - 0(130-140). 



III. SALVAGING MSSM MODELS WITH OVER- ABUNDANT RELICS 

In this section we calculate the degree of relic density dilution needed in MSSM models 
as a function of the temperature at which neutralinos freeze out. This will allow us to 
immediately determine whether a given MSSM model can or cannot be salvaged by a strongly 
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Parameter 


Lower Lim. 


Upper Lim. 


Scan Type 




66 GeV 


20 TeV 


Loot 


|Mi 1 


40 GeV 


20 TeV 


Loe: 


\ \ 


83 GeV 


20 TeV 


Lo£[ 


771, n 


-minfMi Mo a) 


20 TeV 






-3 


3 


Lin 


TUA 


200 GeV 


20 TeV 


Log 


tan/3 


2.5 


60 


Lin 



TABLE L Ranges for the parameter space scan we employ for our fig. d 

first-order electro- weak phase transition, the requirement being that the critical temperature 
of the phase transition be smaller than the freeze-out temperature, and the dilution factor 
be large enough to bring the neutralino relic density at or below the level of the inferred 
average dark matter density in the universe. 

We work in the context of the i?-parity, flavor and CP-conserving MSSM, we enforce that 
the lightest supersymmetric particle be the lightest neutralino, and parameterize the soft 
supersymmetry breaking parameters that enter the relevant particle spectrum for the calcu- 
lation of the neutralino relic abundance with their values at the electro- weak scale (we thus 
do not assume any grand unifled structure for the soft supersymmetry-breaking terms). We 
do not assume any relationship between the gaugino soft supersymmetry breaking masses, 
nor about their relative signs, but we assume a common mass scale for all sfermions. Also, 
for simplicity we set to zero all trilinear scalar coupling with the exception of third gen- 
eration sfermions. The gluino mass (which does not enter into the calculation of the relic 
abundance in the DarkSUSY code) was set to M3 = 3M2, which approximately follows the 
usual supergravity relation. Details of our scan procedure, including lower and upper limits 
for the scan as well as whether the sampling was carried out logarithmically or linearly, are 
given in Table [TTIl 

Once a particular MSSM setup is deflned by the random parameter values picked by the 
procedure outlined above, we require that the resulting particle setup be compatible with 
updated versions of the limits from collider searches, rare decays and electro-weak precision 
measurements described in Ref. [34]. We then calculate the thermal relic abundance of the 
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FIG. 1: A scatter plot of MSSM models on the plane defined by the ratio of the models' thermal 
neutralino relic abundance (Tiparticiey' over the universal dark matter density (^^stro), versus the 
neutralino freeze-out temperature Tf^. Red points correspond to bino-like lightest neutralino, while 
blue and green points to wino- and higgsino-like lightest neutralinos, respectively. Models above 
the horizontal line at f^particie/^astro = 1 ^^^e over- abundant, and are ruled out unless a dilution 
mechanism such as the one we discuss in the present study is operative. See the text for definitions 
and details. 

lightest neutralinos, in the context of a standard cosmological setup, with the DarkSUSY 



package [3J]. 



The two parameters we are interested in for the present study are the thermal neutralino 
relic abundance, and the temperature at which the neutralino freezes out: a phase transition 
occurring at temperatures larger than the freeze-out temperature would not affect the relic 
abundance of neutralinos, since the latter would be in thermal equilibrium after the phase 
transition and its number density would re-equilibrate to the other thermal species. We 
define the freeze-out temperature according to the prescription of Ref . [35] , that sets it as 
the temperature where the comoving neutralino number density is a factor 2.5 larger than 
its asymptotic zero-temperature value. 
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Fig. d] shows the minimal dilution factor required to bring the thermal neutralino relic 
abundance below the upper limit set by determinations of the average universal dark matter 
density, f^astro^^ 0.113 [3] as a function of the neutralino freeze-out temperature. Loosely, 
the latter is a factor of 20-25 times smaller than the neutralino mass, which thus increases 
linearly with the x-axis. We indicate with blue dots models where M2 < Mi, /x, i.e. models 
corresponding to a dominant wino component in the lightest neutralino mass eigenstate; 
green dots correspond to models with < Mi, M2 (higgsino-dominated lightest neutralinos) 
and red dots to Mi < M2, /x (bino-dominated lightest neutralinos). 

Notice that the pair annihilation of bino-like neutralinos proceeds through a variety of 
channels, including squark and higgs exchanges, which dramatically depend on the details 
of the spectrum of the relevant particles. This induces a wide scatter in the relic abundance 
versus mass (or freeze-out temperature) of bino-like models, as evident in fig. [H On the 
other hand, for dominant ly wino- and higgsino-like neutralinos, the dominant annihilation 
modes always proceed through gauge boson pairs, mediated by chargino (or by the next-to- 
lightest higgsino in the case of ZZ final state) exchange. In this case, the pair-annihilation 
cross section is fixed by gauge couplings and by the mass of the neutralino/chargino system, 
modulo kinematic threshold effects for neutralino masses near the gauge boson mass (this 
shows up in fig. [1] for Tfo 2 — 5 GeV). This is the reason why the wino and the higgsino-like 
models all fall approximately on a line in the log-log plot of fig. [H 

MSSM models whose relic abundance can be salvaged by a strongly first-order electro- 
weak phase transition therefore feature either higgsino-like lightest neutralinos with masses 
larger than a TeV, or wino-like neutralinos with masses in excess of 2 TeV, or bino-like 
neutralinos in a wide mass range, provided the freeze-out temperature is large enough to 
be above the temperature at which the phase transition occurs. We discuss the possible 
eflFects of entropy dilution in the electro-weak phase transition and its impact on the MSSM 
parameter space in the following section, where we show which regions of fig.[T]can potentially 
be salvaged by the resulting entropy injection. 
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IV. RELIC DENSITY DILUTION AND FIRST-ORDER PHASE TRANSITIONS 



A. Overview of the Field-Theoretical Setup 

We calculate here the dilution of a thermal relic due to a first-order phase transition 
driven by a scalar field (j). Since we are interested in the specific case of the EWPT, as an 
illustrative instance we consider here the eflFective potential of the neutral component of the 
scalar electro-weak sector of the Standard Model (j). The extension to the Standard Model 
we then take into consideration will include an additional set of fermionic and/or bosonic 
degrees of freedom, which might be thought of as being singlets under the Standard Model 
gauge group, but that might transform nontrivially under a hidden sect or gaug e group, and 



2d, 



2l| . These class 



are coupled to the visible sector only through the Higgs sector 18, 1^. 
of models includes so-called Higgs Portal scenarios, where (partly) secluded hidden sectors 
communicate with the Standard Model (the "visible" sector) exclusively via interactions 
in the Higgs sector \^\. Early studies on such models include those listed in Ref. js?]. 
where the Higgs portal was often invoked in the context of identifying a viable particle dark 
matter candidate. The existence of hidden sectors that might be only partially secluded 
is ubiquitous to extension to the Standard Model, such as those based on high-rank GUT 



models, e.g. where the GUT group is 



38|], in gauge-mediated supersymmetry breaking 



setups (for a review, see e.g. [39|]) as well as in the landscape of string compactifications 
[id ]. In all of these setups, a potentially large number of (fermionic or bosonic) degrees 
of freedom, neutral under the Standard Model gauge interactions, might interact with the 
Higgs sector alone. 

At tree-level, the effective potential associated to the field (j) only we consider here is: 



while the zero-temperature one- loop corrections read: 



^i(0) = E± 



(13) 



(14) 



where gi is the number of degrees of freedom of each particle species in the theory coupled 
to 0, mi{(j)) is the particle species mass, and the upper and lower signs correspond to bosons 
and fermions, respectively. We consider only scalar bosons in this paper. Notice that the 
potential for vector bosons would carry a constant of | in place of the | in Eq. ( |T4l ) (see 
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Ref. {3]). Note that both Vb((/>) and Vi{(l)) have stationary points at = and (j) = v. 
These points are respectively local maxima and minima of Vo((/>). Vi{(j)) has a saddle point 
dit (j) = V and either a maximum or a minimum at = 0, depending upon the leading sign. 
We obtain the free energy density by adding the finite-temperature one-loop correction 



mi{(j)) 
T 



+ E ^ - + n.(r)]^/2] , (15) 



bosons 

where /_ and are for the relevant thermal distribution functions for the bosonic and for 
the fermionic contributions, respectively: 



dy y"^ log (^1 T exp" 



(16) 



and where the second summation is for the resummed Daisy diagrams with n^(T) = |/i^T^, 
where hi is the Yukawa coupling. Assuming all particles acquire mass through a Higgs-like 
mechanism in which the mass terms are of the form m^((/)) = hicj) (i.e. neglecting explicit 
mass terms for all the additional degrees of freedom interacting with the field 0), the free 
energy density takes the form 



647r2 



log- 



^ 27r2 



V 

hi4> 
T 



3 
2 

+ 



.2 j2 



+ 2v 

7r(-lTl 
12 



■D 



hi(j) 



(17) 



where D{x^h) = (x^ + — x^. For our purposes, we can ignore a constant vacuum 

energy term. It will also be interesting in our analysis to add a temperature independent 
cubic term VqX^) = o^(|'^^0^ — ^vc/)^) with a <C A to account for example for the possible 
effects, at the level of the cp effective potential, of tree level cubic terms, driven by mixing 
with one or multiple gauge singlet scalar fields (see e.g. Ref. [41]). Notice that the purpose 
of the quadratic term in Vq is to ensure that the vacuum expectation value is maintained at 
(j) = V dit zero temperature. 
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FIG. 2: The integrals of Eq. lfTS\). Bosons correspond to (black line), and fermions correspond 

to (red line). The dashed lines show the resummed Daisy diagram contributions for h —0.2, 

0.4, 0.6, 0.8, and 1.0. 

In the high-temperature (low-0) hmit, Eq. f fTGl ) can be expanded as 

4 2 4 2 

^ , , TT TT 9 TT o X , X 

'-'"'=-45 + 12" -r -32'°^* 



(18) 



(/ + 2)! ^ ^'\27^J 



/+(x) = 1 X H log — 

+ ^ ' 360 24 32 ^a/ 



+4-' Z.(-i)77^(^i-^Jr(/ + ,)y 



(19) 



(see e.g. the seminal work of Ref. [29]; see also [28]) where loga^ = | — 27 + 21og(47r), 
log a/ = I — 27 + 21og7r, 7 is the Euler constant, C is the Riemann zeta function, and F is 
the Gamma function. 

In the low-temperature (high-(/)) limit, the expansions are given by 

I^{x) = -x' J2 ^-^K,{kx), (20) 
k=i 



(again, see e.g. Ref. [29] and [28|]) where K2 is the modified Bessel function of the second 
kind of order 2. 

There are a few important things to note about these expansions. The original integrals 
in Eq. ( fT6l ) are negative monotonically increasing functions with lima^^oo ^±(^) = 0. They 
are similar in shape to upside-down bell-curves (see Fig. [2]). The low-temperature expansions 
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have this same character for any finite number of terms. If we keep only terms up to A: = 8 
(as we do for the numerical study we present in this analysis), the errors at x = are only 
0.01% and 0.05% for and /_, respectively. The errors drop exponentially for x > 0, 
with fractional errors of ~ 10~^° at x = 2. In contrast, the high-temperature expansion 
diverges towards ±oc for a finite number of terms. For example, retaining even up to / = 15 
produces a visible divergence in towards negative infinity at x = 3.5. The benefit of 

the high-temperature expansions is that they contain an explicit cubic term, whereas the 
low-temperature expansions do not. If the phase transition is weakly first-order so that 
terms of order x^ can be ignored, then the cubic term is necessary to avoid a continuous 
phase transition. Therefore, the low-temperature expansions underestimate the strength of 
the transition when the transition is weak. 

As argued in Ref. {2])], higher-order corrections to the potential are subdominant and can 
be safely neglected. In particnlar, Raf. 0. which adopts a setup vary similar to ours, finds 
that two-loop effects, as caluclated e.g. in [42], do not affect significantly the structure of 
the effective potential. We thus neglect them here. 



B. Semi- Analytical Results in the High- Temperature Limit 

If the phase transition is weakly first-order, we may use a temperature-dependent quartic 
potential as an approximation to the free energy: 

^(0, T) = D{T' - T2)02 + (lav - ET)4>' + - ^giT\ (21) 

The constants and Tq come directly from the coefficients of Eq. (fTTl) and Eqs. (11811191) . 

and gi = gb + ^df are the effective number of degrees of freedom of particle species. For this 
analysis we assume that D is positive, which does not necessarily follow from Eq. f fTTI ). 
It is much more simple to analyze this equation if we recast it in the dimensionless form 

:F{ip, t) = D{t^ - l)ip^ + E{x - T)ip' + IX^' - ^giT\ (22) 

where = J-'/Tq, (p = (P/Tq, t = T/Tq, and x = ^av/ETo. The critical temperature (the 
temperature of degenerate minima), is then easy to find analytically: 



'c 



y-i 
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where y = DX/E^. In all situations, the solution corresponding to the positive root is the 
correct physical solution, while the negative root is an unphysical solution resulting from 
the use of an approximate potential. Note that if x > 1 then by Eq. (I23l) x > Tc and the 
second minimum is at 99 < at r = Tc. Since the temperature-dependent terms in Eq. (fTTl) 
are even in (j) while those in Eq. (I2TI) are not, we reject all solutions with < and demand 
that X < 1. 

If the transition temperature is close to the critical temperature, as we expect for a weakly 
first-order phase transition, then Eq. ([3]) gives the correct amount of dilution. The entropy 
density is 5 = —dT/dT s/Tq = —dJ^ /dr. In the high-temperature phase at the critical 
temperature, we have s^/T^ = ^gir^. In the low temperature phase, the minimum is 
LPc = '^E{rc — x)/\ and the entropy difference is As/Tq = Lpl{2Drc — Eipc)- Plugging in 
values, we get 



A^^J^SE^ y-x^{x^ + y-l)y \ y - x,/{x^ + y - l)y \ 

5+ 2gi7r^ y x^ + y jy x^ + y J ' ^ ^ 

which simplifies to 

As 45 8DE^ 

in the limit that ^ 3> x, 1. Note that the cubic term does not have a noticeable effect upon 
the dilution until it becomes quite large, which would violate our original assumptions. 

In passing, we remark that, taking the above analysis at face value, in the Standard 
Model, where the top quark and the W and Z bosons are the only relevant degrees of 
freedom (see Ref. j22]), and = 0.16, = 0.0096, y ^ 300 (dependent upon the Higgs 
mass), and gi = 106.75, one finds a negligible dilution factor, namely: 

fafV 5+ As 45 8DE^ , , 

M =^^1 + — = 1 + - -——^ 1.0001. (26) 

\CLi J s_ 5+ 2gi7r^ 

Interestingly, the semi-analytic setup outlined above also allows us to get an estimate 

of the dilution factor expected in the context of the minimal super symmetric extension of 

the Standard Model with light scalar tops (see e.g. Ref. |23, 24, 26]), which has often been 

considered in the context of electro- weak baryogenesis. In that case, the cubic term E can 

be one order of magnitude larger than in the Standard Model, implying a dilution factor of 

at most 1.01. According to our equation f l25l ) this means that in the MSSM the dilution to 

the relic density of species freezing out prior to the electro- weak phase transition is a small 

effect (at most a few percent). 
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FIG. 3: Dilution factors for boson-fermion models (dotted red and dashed green) and purely bosonic 
models (solid black) as functions of the temperature at the end of reheating (left) and the number 
of fermion and boson degrees of freedom (right). For the purely bosonic models, the highest line 
corresponds to the highest Yukawa coupling for which the point at (j) — Q is a minimum at zero 
temperature. The three lower lines have 0.7, 0.8 and 0.9 times this value. For boson-fermion 
models, the four dotted red lines correspond to Yukawa couplings identical to those in the bosonic 
cases, while the dashed green lines corresond to Yukawa couplings of 0.5, 0.75, and 1.0. 

C. Strongly First-order Phase Transitions: Numerical Solutions 



If the phase transition is strongly first-order, Eq. f l2TI ) does not hold and we must instead 
resort to the full expression for the free energy density with the low-temperature expansion 
of Eq. (fT6l) . One must compute the dilution factor fully numerically. 

In our models, we proceed with the following steps. First, we obtain a function of the 
minimum of the true vacuum 00 (r) up to the critical temperature by numerically solving 
the differential equation 

* = _ (^\ I (^2L\ (27) 

with the initial condition (/)o(T = 0) = v. We then calculate the transition temperature by 
calculating the Euclidean action at many temperatures and searching for S-^/T ^ (9( ISO- 
MO). To account for a potentially wide range of transition temperatures and relativistic 
degrees of freedom, we use Ss{T^)/T^ = 170 — 51og(T/l GeV) — 21og(5'). We then calculate 
the dilution factor directly from Eq. ([9j) using s = —dT/dT. 

For definiteness, we consider two classes of models: one with only additional scalar bosonic 
degrees of freedom, and one with an equal number of additional bosonic and fermionic 
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FIG. 4: Reheating temperature as a function of fermion and boson degrees of freedom. Lines cor- 
respond to the same models as in Fig. O with higher Yukawa couplings resulting in lower reheating 
temperatures. 



contributions. An explicit realization o 
secto. scala. S, d.seu.ed in Re. U 



le former case is given by the multiple hidden 



2l| . where the coupling between the hidden and 



the visible sector occurs only through renormalizable terms in the potential proportional to 
H^HSf. The latter, instead, alludes to a supersymmetric particle content for the additional 
degrees of freedom; one possible example is an extension to multiple extra generations of 
fermion-sfermions along the lines of the analysis of Ref . [43] , that discusses the impact on 
the EWPT of a supersymmetric model with four chiral matter generations. 

For the ease of analysis, all massive particles in a given model have the same Yukawa 
coupling so that we effectively have only one (or two) massive particle species with a poten- 
tially large number of degrees of freedom. In each case, the degrees of freedom per fermion 
and/or boson range from 10 to 10,000. We additionally add massless bosons (i.e., h = 0) 
with gi = 100 degrees of freedom, which approximates the number of relativistic degrees of 
freedom in the standard model. That is, we have a variable number of degrees of freedom for 
massive particles, but keep the number of massless particles constant. The latter just adds 
the ^-independent term —7r'^giT^/90 to the free energy of Eq. f fTTI ). and gives an additional 
suppression to the factor ( of Eq. f fT2l ). Note that the addition of massless particles does 
not change the phase transition temperature or dynamics. It does, however, add entropy to 
both the high and low-temperature phases, which decreases the overall dilution factor via 
Eq. [9l In the boson-only models, large couplings turn the point at = into a maximum 
at zero temperature via the quadratic term in Eq. ( iMl ) . We limit ourselves to the maximum 
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couplings for which this is not the case. In aU cases, we set the zero-temperature Higgs mass 
to be rrih 150 GeV. 

Figures [3] and H] display our results. The black lines correspond to models without any 
fermion contribution and variable Yukawa couplings. The models in the highest of these 
lines have the maximal couplings for which the point at = 0, T = is a maximum, while 
the lower lines have 0.9, 0.8, and 0.7 times the maximum couplings. The dotted red and 
dashed green lines correspond to models with equal fermion and boson contributions. The 
dotted red lines have the same couplings as the solid black lines, while the dashed green 
lines have fixed couplings of 0.5, 0.75, and 1.0. For boson- fermion models with reasonably 
large Yukawa couplings and many degrees of freedom, one can achieve dilution factors on the 
order of 100. Purely bosonic models, on the other hand, can only reach dilution factors on 
the order of 10. The largest aspect contributing to this discrepancy is the range of Yukawa 
couplings available to the diflFerent cases. At g = 10,000, the maximum coupling for bosons 
is only h^, = 0.23, while the highest plotted dilution factor for a boson- fermion model has 
hi) = hf = 1.0. When purely bosonic and boson- fermion models have identical couplings, 
the bosonic models tend to have slightly higher dilution factors. 

One can easily explain two of the important qualitative features of Figs. [3] and [H by 
examining the free energy of Eq. (fTTI) . First, the temperature of the phase transition tends 
to decrease with increasing degrees of freedom. The temperature-dependent term in the free 
energy has a coefficient of gT^. When this term is large, the second minimum of the potential 
disappears and the origin becomes the true vacuum. Thus, we expect the temperature of 
the transition to scale as T oc g~^^^^ as seen in Fig. [H Second, the dilution factor tends to 
increase with increasing Yukawa couplings and decreasing temperature scales. If either h 
is large or T is small, then the temperature-dependent term is nearly zero except for small 
(/), excluding resummed Daisy terms (see Fig. [2] for the behavior of I±{x)). Therefore, the 
entropy s = —dT/dT is much smaller in the true vacuum than it is at the origin. Including 
the Daisy terms, the negative slope of — ^D{x^h) tends to decrease the entropy at 

high (/). If the transition reaches a phase coexistence stage, then the small entropy leads 
directly to a high dilution factor via Eq. ([3]). This effect dominates when h(j)/T ^ hv/T 5^ 
which is the case for the boson- fermion models in Fig. [3] with h > 0.5 (dashed green lines). 
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FIG. 5: Dilution factors for models with cubic terms (red lines) added to the boson-fermion models 
in Fig. with Yukawa couplings of h — 0.5 (black line). The red lines have cubic strengths of 
a/X= 0.05, 0.15, and 0.50. 

D. Models with Tree-Level Cubic Terms 

We also examine what happens when we add a tree-level cubic term V'o(^) = oi{\v^(lP' — 
^v(l)^) to the free energy. Fig. [5] shows cubic terms added to boson-fermion models with 
Yukawa couplings fixed at /i/ = /15 = 0.5. The models are otherwise exactly the same as 
those in the previous section: gi = 100, rrth = 150 GeV, and the boson/fermion degrees of 
freedom gf = g^, range from 10 (corresponding to the lower-right portion of Fig. [5]) to 10,000. 
The black line has no added cubic term, while the red lines have cubic strengths of a/A = 
0.05, 0.15, and 0.50. Even large terms with a = |A do not seem to significantly impact 
the dilution factor. The effects of cubic terms upon purely bosonic models and models with 
smaller Yukawa couplings are qualitatively similar. 

V. SUPERSYMMETRIC DARK MATTER AND A STRONGLY FIRST-ORDER 
PHASE TRANSITION 

We determine here which MSSM models with neutralino dark matter might not overpro- 
duce dark matter given a strongly first-order phase transition. Figure [6] shows an overlay of 
the dilution from diflFerent models of phase transitions on top of our scan of the dark matter 
abundance in MSSM models. Red, green and blue points correspond to bino, higgsino, and 
wino-like models, while the solid black, dotted red, and dashed green lines correspond to 
phase transition models with only bosons at their maximal Yukawa couplings, bosons and 
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FIG. 6: An overlay of dilution factors for boson and boson-fermion models on top of a scatter plot 
of dark matter abundance in different MSSM models. The points are the same as those in Fig. d 
The green, black, and red lines correspond to the highest green, black, and red lines of Fig. The 
two boxes show which MSSM models are viable for dilution models with hf — hij — 1.0 and Qh^Qf < 
10,000 and 1,000. 

fermions at the same Yukawa couplings, and bosons and fermions with fixed Yukawa cou- 
phngs of /i/ = /15 = 1.0. These are the same as the highest of the black, red and green lines 
in Fig. [31 

The freeze-out temperature of the dark matter must be larger than the final reheating 
temperature of the phase transition if the transition is to have an effect upon dark matter 
abundance. Therefore, any point that lies to the right of any of the three lines can be 
diluted to the cosmic relic abundance by a sufficiently strong phase transition. If a point 
lies between the lines, then one can easily pick a phase transition model that produces the 
correct amount of dilution. Note that only MSSM models with an overabundance of dark 
matter are problematic — if a given phase transition model over-dilutes a given MSSM model 
of dark matter, then the two models can still be mutually viable if the neutralino dark matter 
does not constitute the entire universal dark matter abundance. 

The two boxes in Fig. [6] show which MSSM models are viable given the most optimistic 
phase transition models (that is, boson-fermion models with /ij = /15 = 1.0) with degrees of 
freedom gf = gt = 1,000 (inner box) and 10,000 (outer box). With 1,000 d.o.f., almost all of 
the over-abundant higgsino and wino-like models are viable. With 10,000 d.o.f., a handful of 
bino-like models become viable as well. However, there are no viable models with freeze-out 
temperatures less than ^ 20 GeV or overabundances greater than a factor of ^ 200 unless 
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they are accompanied by phase transitions with very large (> 10,000) particle degrees of 
freedom. This includes most of the bino-like models and a small subset of the low-mass 
higgsino-like models. 

VI. SUMMARY AND CONCLUSIONS 

A strongly first-order electro- weak phase transition, warranted in the context of scenarios 
for the production of the observed baryon asymmetry at the electro-weak scale, can lead to 
a significant dilution of the thermal relic abundance of the lightest supersymmetric particle, 
if the latter is stable. In this paper we studied the impact of toy models for the electro- weak 
sector on the relic abundance of MSSM neutralinos. Specifically, we considered the effect on 
the electro- weak sector of adding a large number of bosonic or bosonic-plus-fermionic degrees 
of freedom, which affects both the one-loop temperature-independent and the temperature- 
dependent contribution to the effective potential of the Higgs field. We noted that for 
neutralino freeze-out temperatures between 20-40 GeV, corresponding to neutralino masses 
at or above 400 GeV, the dilution due to the entropy injected in the first-order electro- weak 
phase transition can be as large as 10-100 if we postulate a large number of extra degrees of 

freedom (10^ 10^). Numerous MSSM models exist that could be viable if such a dilution 

effect is in fact operative. For models with a more modest number of additional d.o.f. (say, 
of the same order of the Standard Model d.o.f.), the dilution can still be on the order of 2. 
We find that a cubic term in the effective potential has a comparatively small effect upon 
the relic abundance of species. Also, we showed that the dilution expected in the context 
of the MSSM without additional degrees of freedom coupled to the Higgs sector is a very 
small effect, at most at the few-percent level. 

Should future data from the LHC point towards a neutralino relic abundance larger than 
the cosmological dark matter density, the origin of the mis-match might lie in the same 
electro-weak scale physics that the LHC itself will concurrently explore. Interestingly, this 
can be profoundly intertwined with the question of the origin of the matter-antimatter 
asymmetry in the early universe, opening up the possibility that electro-weak physics lie at 
the core of both dark and visible matter, and that the LHC will soon shed light and perhaps 
unveil this scenario. 
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